Hybridized solid-state qubit in the charge-flux regime 
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Most superconducting qubits operate in a regime dominated by either the electrical charge or the 
magnetic flux. Here we study an intermediate case: a hybridized charge-flux qubit with a third 
Josephson junction (JJ) added into the SQUID loop of the Cooper-pair box. This additional JJ 
allows the optimal design of a low-decoherence qubit. Both charge and flux 1// noises are considered. 
Moreover, we show that an efficient quantum measurement of either the current or the charge can 
be achieved by using different area sizes for the third JJ. 



I. INTRODUCTION 

Solid-state qubits based on Josephson-junction (JJ) 
circuits have attracted considerable attention in recent 
years and different kinds of Josephson qubits are being 
explored by taking advantage of the charge and phase 
(flux) degrees of freedom. Experimentally, quantum os- 
cillations were observed in chargcji phase,^ and flux 
qubits."^ A Josephson qubit in the intermediate regime 
between charge and flux also exhibited quantum oscilla- 
tiona^ and showed a high quality factor corresponding to 
a decoherence time of about 0.5 fis. Because quantum in- 
formation processing requires states to evolve coherently 
in a sufficiently long time, it is thus crucial to obtain 
qubits with very low decoherence. 

Here we study a new type of Josephson qubit, some- 
what similar to that in Ref. 0, in which a third JJ 
is added into the SQUID loop of the Cooper-pair-box 
(CPB) qubit. ^ In Ref. 4 this third JJ is connected to a 
current source and only used for measuring the quantum 
states of the CPB qubit. Moreover, because this JJ is so 
large, the quantum states of the CPB qubit are only very 
slightly modified by it. Actually, without a bias current, 
the large third JJ can be approximated by a harmonic 
oscillator and the whole system can thus be considered 
as a CPB qubit coupled to the oscillator. This is very 
similar to a CPB qubit in a cavity^M. Very recently, the 
coherent dynamics of a flux qubit coupled to a harmonic 
oscillator has been studiedji where the large-size SQUID 
connected to the qubit plays the role of the harmonic 
oscillator. 

In our present work, the CPB qubit is working in the 
charge-flux regime, as in Ref. 0, but now the third JJ is 
not necessarily large and more importantly it is not used 
just as a measuring component. This additional degree 
of freedom in designing the charge-flux qubit allows us 
to optimize the qubit by changing the size of the third 
JJ. Indeed, here we show that the charge-flux qubit is 
gradually hybridized (in the quantum mechanical sense) 
with the third JJ when the area size of this additional 
JJ decreases. More importantly, we find that the qubit 



can be optimized to have the lowest decoherence at a 
suitable size of the third JJ. Furthermore, we show that 
efficient quantum measurements of either the current or 
the charge can be implemented by just choosing different 
sizes for the third J J. 

The paper is organized as follows. In Sec. II, we present 
the model Hamiltonian for the hybridized charge-flux 
qubit and study its properties. Section III shows the 
energy spectra of the qubit for different sizes of the third 
junction added into the SQUID loop of the CPB. To an- 
alyze the effects of different kinds of noises on the qubit, 
we employ the boson bath model in which a noise is de- 
scribed by a collection of spectrally distributed harmonic 
oscillators. The characteristic times for relaxation, deco- 
herence, and leakage of the qubit states are calculated in 
Sec. IV. We optimize the qubit to have the lowest deco- 
herence by choosing a suitable size for the third junction. 
Section V is devoted to quantum measurement. We pro- 
pose two readout schemes to efficiently discriminate qubit 
states by taking advantage of the charge and flux degrees 
of freedom. Finally, conclusions are presented in Sec. VI. 



II. THE MODEL 

The hybridized charge-flux qubit is shown in Fig. 1(a). 
The third JJ, i.e., the left one, is added into the SQUID 
loop of the CPB in which an island (denoted by a black 
dot) is connected by two JJ and coupled to a gate voltage 



by a capacitance Cg. When Ci 
nian of the system is given by 



Ca = Co , the Hamilto- 



H ^ E,p{N - ngf + Ei{N3 + -ngf + U, (1) 



with 



Here 



U ^J2E,n{l-cos<j>,). 



Ecp = 2Ec, Ec — 

ZLj2 



(2) 
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El = 



8C2E12 
C2 + 2C3 



(3) 



The phase drops through the three junctions are con- 
trained by 



where 



bs + 27r/e = 0, 



f 



(4) 



(5) 



is the reduced magnetic flux in the qubit loop (in units 
of the flux quantum $0 — h/2e). The operator 



N = -i 
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corresponds to the number of Cooper pairs on the island, 
and 



d 




a = 20 




no left junction 
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corresponds to the number of Cooper pairs tunneling 
through the left J J. Here we consider the simpler case 
with Eji = Ej2 ^ Ej, Ci ^ C2 ^ C, E,,3 = aEj, and 
C3 = (3C. In this case, the periodic potential ^a) is 

U ::=-Bj[(2 + Q)-2cos0cos(7r/e + i(/;>3)-acos(/)3], (6) 

and the condition Ci + Cg — C2 can be approximately 
achieved because Cg <C Ci , C2 . 

Assuming that the eigenstate of Hamiltonian has 
the form as follows 



^{(t>, 03), 



(7) 



one can cast the equation for eigenvalues and eigenfunc- 
tions 



i)^E<^{ 



(8) 



to a standard Schrodinger equation with a periodic po- 
tential: 



where 



FoV'(</',03) ^i^V-C-^^s), (9) 



Ho = E.pN'' + EpNl + U{^, 03), (10) 



with U{(j),(j)3) given by Eq. 

Similar to a flux qubit (see, e.g., Refs. H and0), the 
reduced Hamiltonian Hq is just like that for a particle in 
a two-dimensional periodic potential, so the solution of 
Eq. © has the Bloch -wave form 



(11) 



FIG. 1: (Color online) (a) Schematic diagram of the hy- 
bridized charge-flux qubit, which consists of three JJs in a 
superconducting loop (pierced by an external magnetic flux 
"I>e) and a superconducting island (denoted by a black dot) 
coupled to a gate voltage via a capacitance Cg. The Joseph- 
son energies and capacitances of the JJs are Eji = Ej2 = Ej, 
Ci = C2 = C, Ejz = Q-Ej, and C3 = /3C. Here we consider 
the charge-flux regime when Ej ^ Ec = e^/2C2. Unless 
explicitly stated otherwise, q = /3 is chosen throughout the 
paper. Energy levels of the charge-flux qubit versus (b) rig at 
/e = 0, and versus (c) /e at Ug — 0.5, where a — f3 — 20. (d) 
Energy levels of the CPB qubit versus rig at /e = 0, without 
the left JJ. 



where K = (fc, /C3). The constraint 

{k,k3) = i-ng,^ng) 
on the wave vectors gives rise to 

«-(?!), 03) =Uk(0, 03), 



(12) 



(13) 



which ensures that 5'((/), ^3) is periodic in the phases 
and 03. 

Moreover, the Hamiltonian |^ can be rewritten as 



H — Hep + Hi + Hi, 



(14) 



where 



H,p = E,piN - ngf + 2E.,[\ - cos0cos(^/e)] (15) 

is the Hamiltonian of a CPB qubit, i.e., the qubit with 
the left J J absent in Fig. 1(a), and 

Hi = Ei{Ni + ^ng f + aEj{l - cos 03) (16) 



3 



is the effective Hamiltonian of the left J J. The interaction 
Hamiltonian 

Hi = 2Ej cos cf) [cos(7r/e) - cos(7r/e - '^(p^)] (17) 

represents the couphng between the CPB qubit and the 
left J J. 

For a large left JJ, the phase drop 03 is small, so the 
left J J can be approximated as a harmonic oscillator with 
frequency 



n = -{KEjE,f'\ 



where 



a 



(18) 



(19) 



Also, the interaction Hamiltonian can be approximated 
by 



Hi 



A3 sin(7r/e) - cos(7r/e)]£',/ cos 0, (20) 



with 



4& 



-,1/4 



a{l + f3)Ej 



(a + a^). 



(21) 



where a {a^) is the operator for annihilating (creating) 
a boson. Because a and /3 are large for a large-area left 
JJ, it is clear that when /e ^ 0, Hi is dominated by a 
weak one-boson process, while a much weaker two-boson 
process is involved in Hi for /e = 0. 



III. ENERGY SPECTRUM 

Below we show the hybridizing effects of the left JJ 
on the energy spectrum of the qubit in the charge-flux 
regime with Ej = Ec- The energy levels for /e = and 
Ug = 0.5 are given in Figs. 1(b) and 1(c), where a large 
left JJ with a = /3 = 20 is chosen. In contrast to the 
energy levels of the CPB qubit [cf. Fig. 1(d)], there exist 
additional levels due to the left JJ. However, because the 
left JJ is now large (i.e., Ej^ — 20Ej), the interaction 
between this JJ and the CPB qubit is small. Therefore, 
the energy levels of the CPB qubit are slightly modified 
by these additional levels, especially for the two lowest 
levels used for the qubit. 

Figures 2(a) and 2(b) display the energy levels for 
fe — and Ug — 0.5 and a much smaller Ejs, since 
now a = 3. The levels of the left JJ now hybridize with 
those above the two lowest levels, but the two lowest 
levels are still barely modified [comparing Fig. 2(a) with 
Fig. 1(d)]. This means that, as far as the two lowest 
states are concerned, the left JJ with a = 3 can still be 
regarded as a large J J. When the left J J becomes even 
smaller (e.g., a = 0.3), Hi becomes larger and the energy 
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FIG. 2: (Color online) Energy levels of the charge-flux qubit 
versus rig at /e = for (a) a = 3 and (c) 0.3, and versus /e 
at Ug = 0.5 for (b) a = 3 and (d) 0.3. 



levels of both the left JJ and the CPB qubit become 
heavily hybridized [see Figs. 2(c) and 2(d)]; one can see 
that the energy levels in Fig. 2(c) look different from 
those in Fig. 1(b), but the two lowest levels can also be 
used for a qubit. 



IV. STATE COHERENCE AND QUBIT 
OPTIMIZATION 



Realistic qubit circuits will experience fluctuations 
from both charge and magnetic flux. These noises will af- 
fect the coherence of the qubit states in the subspace with 
basis states |0) and |1), corresponding to the two lowest 
levels. To characterize the qubit-state coherence, the re- 
laxation time Ti and decoherence time T2 are usedtifi 



1 

1 



1 



with 



- = |(0|A|0)-(l|A|l)|25(c.)|^^o. 



(22) 



(23) 



Here A is an operator characterizing the coupling be- 
tween the qubit and the environment, and S{u)) is the 
power spectrum of the noise. Moreover, because there 



4 



are other levels above the lowest two, leakages from the 
qubit-state subspace to these outside levels can occur. 
Therefore, two additional timesti^ 



1 



iY.\{n\A\k)\'SM, 

n 

k = 0,l, n = 2,3,---, 



(24) 



are needed to characterize the noise-induced transitions 
from the two lowest levels to the ones above. 

These results are based on the boson bath model in 
which the noise is described by a collection of harmonic 
oscillators with a spectral distribution. When Eq. 1)22(1 is 
applied to a 1// noise (see, e.g., Ref . . the very low 
frequencies are cut off for the power spectrum S{uj). This 
cutoff low-frequency part corresponds to the limit of very 
slow processes. For instance, for the 1// charge noise, 
this can correspond to the extremely slow switchings of 
the trapped charges. If these fluctuating processes are 
much slower than the decoherence time T2 of the qubit, 
they remain approximately static during the quantum 
operation and yield negligible dephasing. 
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A. Johnson-Nyquist noises 

For Johnson-Nyquist noises, such as the fluctuations 
from gate voltage and external magnetic flux, the opera- 
tors A are given, respectively, by 



Av - 



E,„N - E1N3 



(25) 



and 



A<s> = cos{(f>) sin(7r/e -f -^a)- 



The power spectrum is given by 



S{uj) = J{uj) coth 



tlLU 

2kBT 



(26) 



(27) 



where J{uj) is the bath spectral density. 

For gate-voltage fluctuations characterized by an 
impedance Z{uj), the bath spectral density is 



-Kg 



where 



1 



(1 + 2ay 



£1 
C2 



(28) 



(29) 



and Rq = h/e^ ~ 25.8 kfi is the quantum resistance. 
Here we choose Cg = O.OIC2, and consider the typical 



FIG. 3: (Color online) Relaxation (Ti), decoherence (T2), and 
leakage {Tlo and Tli) times versus Ug at (a) a — 20, (b) 3, 
and (c) 0.3 for gate- voltage noise, and at (d) a = 3 and (e) 0.3 
for flux noise, where fe = 0. In (d), T\ and T2 are not shown 
because they are 5 orders of magnitude larger than Tlo and 
Tli- Here the temperature is chosen to be T = 30 mK. 



Ohmic case of Z{lo) — Ry — 50 fi. The external mag- 
netic flux in the qubit loop is produced by a coil of induc- 
tance L and resistance Rl- The bath spectral density of 
the external magnetic flux fluctuations is 



- 2 y-R-L) [l + {.LlR,n 



where 



MI, 
$0 



(30) 



(31) 



with Ic = 2t:Ej /^Q, and M is the mutual inductance 
between the qubit loop and the coil. Here we choose 
Ej/h = 20 GHz, Rl = 100 f], L = 30 pH, and M = 5 
pH. These parameters correspond to realistic circuits. 

Figures 3(a)-3(c) show the four characteristic times at 
a = 20, 3, and 0.3 for the gate-voltage noise. These four 
times are almost of the same order of magnitude for dif- 
ferent sizes of the left J J; especially Ti (the minimum of 
them) and T2 do not change much. This implies that 
the gate- voltage noise is mainly determined by the ratio 
of Ej/Ec and less sensitive to the variation of the left 
JJ. The observation that Tlo and Tli have almost the 
same order of magnitude as Ti and T2 also means that 
in this case the leakages produce equivalently important 
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effects on the qubit states as the relaxation and deco- 
herence in the qubit-state subspace. However, for the 
noise due to external flux fluctuations, the leakages dom- 
inate over the relaxation and decoherence [cf. Figs. 3(d) 
and 3(e)]. Moreover, when the external magnetic flux is 
around zero, the effects of the external flux noise are sen- 
sitive to the variation of the left JJ. For instance, when 
the left J J decreases in size to a = 3, the leakage times 
in the ffux-noise case are comparable to the relaxation 
time Ti in the case of gate- voltage noise [comparing Fig. 
3(d) with Figs. 3(a)-3(c)], and the qubit-state leakages 
become more serious with a decreasing further [see Fig. 
3(e)]. 

The above numerical results for the Johnson-Nyquist 
noises show that the decoherence time T2 ^ 30 fis or 
more, much longer than the experimental value T2 ~ 
0.5 in Ref. 4. This indicates that they could not 
be the major sources of decoherence in the charge-flux 
qubit. Instead, because the l/f noise may be the main 
source of decoherence, we further study its effects on the 
charge-flux qubit. 



B. l/f charge and flux noises 

There have been numerous attempts to model l/f 
noise; including using a collection of independent bistable 
fluctuators with a given distribution of flipping ratesi^ii^ 
or by interacting two-level classical fluctuators Alter- 
natively, one can also model it using a boson bath with 
a 1// spectral density>i^ For the charge-flux qubit con- 
sidered here, there can be two independent l/f charge 
noises related with the background charge fluctuations of 
the CPB and the left J J; the leakage rates l/T^k as well 
as the relaxation and decoherence rates 1/Ti [i = 1,2) 
are the sum of their respective contributions. These two 
charge noises can be characterized by the power spectra 



2E. 



cp 



a, 



he 
he 

with the corresponding operators A being 

d 



g 

LU 



(32) 



Ag,cp - Ig^, 

■ 9 

d(l)3 



(33) 



Here, for simplicity, a = (3, and aq is chosen here to be 
the same for these two charge noises. In Ref. 0, > 
a because a current source is connected in paralell to 
the left JJ; this decreases Ei = 4Ecp/{l + f3), and the 
dephasing due to the l/f charge noise of the left J J is 
weaker than that oi (3 — a. Also, we can define a power 
spectrum for the l/f flux noise: 



/2ttEj 



V ^^0 



QfcJ) 



(34) 
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FIG. 4: (Color online) Relaxation (Ti), decoherence (T2), and 
leakage (Tlo and Tli) times versus a at {ug, /e) = (0.5, 0) in 
the presence of l/f charge noise, where Ej/h = 20 GHz is 
chosen. Inset: Ti and T2 are replotted with a scaled linearly. 



The corresponding operator A is 



A^ = cos(0) sin(7r/e + -(/>3) 



(35) 



which is identical to Eq. (|26|1 . 

Figure 4 shows the four characteristic times at the de- 
generacy point [ug, fe) — (0.5,0) in the presence of the 
l/f charge noise. We choose aq = (0.7 x lO^'^e)^ for 
the power spectrum, which is very close to the value 
used for fitting the experimental data of the l/f noise 
in the charge qubit The cutoff frequency is chosen 
to be LOc/^-n = 60 Hz, corresponding to a time scale 
~ 2 X 10** /iS, much slower than the experimentally mea- 
sured decoherence time 0.5 /is of the charge-flux qubit >i 
To compare the effects of both charge and flux noises, 
we take the same cutoff frequency for the l/f flux noise. 
Moreover, we use q;$ = 3 x 10~^'^^q for the flux-noise 
power spectrum, which is the experimentally determined 
value of the flux qubitii In Fig. 4, the obtained leakage 
times Tlo and Tli are longer than Ti and T2. This means 
that the leakage is not significant for the l/f charge noise, 
even though the two lowest levels for the qubit are not 
very separated from the higher levels (cf. Figs.l and 2). 

We also calculated the four characteristic times for the 
l/f fiux noise at {ng,fe) = (0.5,0) and found that they 
are much longer than the corresponding characteristic 
times for the l/f charge noise. Thus, we conclude that 
the l/f flux noise plays the least dominant role at the 
degeneracy point for the qubit in the charge-fiux regime 
with Ej — Ec- Moreover, for both l/f charge and flux 
noises, we found that, in the vicinity of the degeneracy 
point (rig, /e) = (0.5, 0), Ti, Tlq, and Tli depend weakly 
on rig and fe- The decoherence time T2 also depends 
weakly on fe (ug) for the l/f charge (flux) noise, but very 
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FIG. 5: (Color online) Circulating currents in the qubit loop 
versus fe for eigenstates \i), (a) i = and (b) 1, where Ug — 
0.5. The number of Cooper pairs on the island versus Ug at 
eigenstates \i) for (c) a = 20, (d) 3, and (e) 0.3. 



strongly on Ug (/e); slightly away from the degeneracy 
point along Ug {fe), the decoherence time T2 decreases 
several orders of magnitudes. 

For clarity, the relaxation and decoherence times Ti 
and T2 are replotted in the inset of Fig. 4 for the 1// 
charge noise. At a — 20, T2 « 0.5 ^s, and Ti « 1.6 ^s. 
Note that this agrees with the experimental results'* of 
the charge-flux qubit with a large left JJ. Also, it can 
be seen that the relaxation time remains at Ti 1.5 fis 
until a 1.3, while the decoherence time T2 first in- 
creases with decreasing a, then remains at r2 3 fis 
(the longest decoherence time) when 1.3 ^ a < 16, and 
finally falls down for a < 1.3. Therefore, one can opti- 
mize the charge-flux qubit in the region 1.3 ^ a 16, so 
that the qubit has the lowest decoherence. 



V. EFFICIENT QUANTUM MEASUREMENT 

Finally, we focus on how to raise the readout efficiency. 
Figures 5(a) and 5(b) display the circulating current / 
in the qubit loop at eigenstates |0) and |1). It is clear 
that the currents {i\I/Ic\i), for a = 20 and 3, are close 
to each other. This further indicates that the left JJ 
with a = 3 still behaves like a large JJ barely affecting 
the CPB qubit. Here we consider the readout scheme 



in Ref. where a current pulse is applied to the qubit 
circuit via a current source connected in parallel with the 
left J J. This gives rise to an effective capacitance C3 with 
a larger value of (3. Thus, the effect of the left J J on the 
CPB qubit is further weakened because the interaction 
Hamiltonian Hj decreases when increasing (3. 

For a single left J J without the right CPB in Fig. 1(a), 
when biased by a current pulse, it switches at 



(36) 



with a narrow switching-probability distribution, from 
the zero- voltage state to the dissipative nonzero- voltage 
state. Here Ic3{o() = ale is the critical current of the 
left J J. However, when the current pulse is biased to the 
qubit circuit, i.e., the left JJ plus the right CPB [see 
Fig. 1(a)], the left J J switches at 



(37) 



with probabilities pi (si) {i — 0, 1) which depend on (see, 
e.g., Ref. 113) 



ic3 



(38) 



In Ref. 0, a = 20, and the switching-probability differ- 
ence is found to be as small as 



because 



bo -Pi I -0.1 



|(0|///.|0)-(l|///e|l)| 



(39) 



(40) 



is small for a = 20. However, when the left JJ becomes 
smaller, to a = 3, 

|(0|///,|0)-(1|///,|1)| 

remains nearly unchanged, but |so — sij is enlarged about 
seven times. This greatly increases \po — pi\ and thus 
efficiently discriminates the states |0) and |1). 

In Figs. 5(c)-5(e), we show the number of Cooper pairs 
on the island, {i\N\i), at eigenstates \i), z = 0, 1. For a 
given i, (i|A''|i) at a = 20 and 3 are similar to each other 
but much different from that at a = 0.3. For instance, 
when Ug = 0.34 (indicated by a vertical black dotted 
line), the number difference 

AiV ^ |(0|iV|0)-(l|7V|l)| 

is AiV ^ 0.37 for both a = 20 and 3, but increases to 
AiV ^ 0.67 when a = 0.3. Therefore, the readout ef- 
ficiency for discriminating the states |0) and |1) can be 
much increased at a = 0.3, when a single-electron transis- 
tor— is capacitively connected to the island and used for 
measuring the quantum states. Also, one can effectively 
couple two CPB qubits with a — 0.3 by taking advantage 
of the charge degree of freedom, such as connecting the 
islands in the two qubits via a mutual capacitance (see, 
e.g., Ref. 0). This capacitive coupling can be used to 
reduce decoherence in a logical qubit composed of two 
CPB qubitsiSa 
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VI. CONCLUSION 

In conclusion, we have studied a hybridized charge-flux 
qubit in which an additional JJ is added into the SQUID 
loop of the CPB. The goal is to find a low-decoherence 
superconducting qubit. This is one of the most important 
open issues since quantum computing is possible only 
if a qubit with long enough decoherence time becomes 
available. 

Currently, 1// noise is believed to be the main source 
of decoherence in a superconducting qubit. Here we con- 
sider the effects of both charge and flux 1// noises on the 
hybridized charge-flux qubit. We find that the qubit is 
optimized in the region 1.3 ^ a ^ 16, so that the qubit 
has the lowest decoherence. These results indicate how to 
optimize a qubit that is expected to have a longer deco- 
herence time. Moreover, we find that the readout scheme 
via measuring currents, like that in Ref. 0, can also be 
optimized, so that the efficiency for discriminating qubit 
states is much increased. Furthermore, we show that 
an efficient readout scheme by measuring charges can be 
achieved as well. 

Note that our studies on the 1// noise use the har- 
monic bath model with a l/w spectral density. This is 
valid when the 1// noise is not dominated by a few fluc- 



tuators strongly coupled to the qubit. Here, due to the 
lack of available data for a charge- flux qubit, the numer- 
ical values of aq and a$ in the power spectra of 1// 
charge and flux noises are chosen from the experimental 
data of the charge and flux qubits. Also, for simplicity, 
the same frequency cutoff is used for both charge and flux 
1/f noises. In some cases, this might considerably devi- 
ate from the realistic samples. Thus, more experimental 
data are needed for giving a quantitative comparison with 
realistic samples. 
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